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Abstract 

On-site measurements of water salinity (which can be directly evaluated from the electrical conductivity) in deep- 
sea sediments is technically the primary source of indirect information on the capacity of the marine deposits of 
methane hydrates. We show the relation between the salinity (chlorinity) profile and the hydrate volume in pores to be 
significantly affected by non-Fickian contributions to the diffusion flux — the thermal diffusion and the gravitational 
segregation — which have been previously ignored in the literature on the subject and the analysis of surveys data. We 
provide amended relations and utilize them for an analysis of field measurements for a real hydrate deposit. 
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1. INTRODUCTION 

Since being discovered methane hydrates have attracted significant attention. First estimates of their amount on 
the Earth, especially in the marine sediments, and their importance were extremely exaggerated. Today's assessment 
of their amount and role are more moderate and well underpinned by field data and results of numerical modelling 
and, thus, may be treated as realistic. Even with this "moderate" evaluation, research on natural hydrates in marine 
sediments is considered to be important. 

In particular, methane hydrates present a potential hazard under anthropogenic climate change. The sensitivity 
of hydrate stability to changes in local pressure-temperature conditions and their existence beneath relatively shallow 
marine environments, mean that submarine hydrates are vulnerable to changes in bottom water conditions (e.g. warm- 
ing). The potential climate impact of methane release following dissociation of hydrate in the past has been compared 
to climate feedbacks associated with the terrestrial biosphere and identified as a possible trigger of abrupt climate 
change (ej*;, IB 01)- The role of hydrate disassociation as a trigger for submarine landslides has also been investi- 
gated Hl|3[|4j], with reports of known hydrate occurrences that coincide with slumping and submarine landslides bein 
common ||J]7 It is therefore imperative to improve our understanding of the global hydrate inventory. Studies J^, 
highlight that this improvement requires also certain revision of the physical and mathematical models of the marine 
deposits of methane hydrates employed in the literature Mam. 
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Figure 1: (Color online) Sketch of the marine sediments hosting a hydrate deposit with a free gas zone (or "bubble horizon") beneath the zone of 
the thermodynamic stability of methane hydrate. The boundary between the hydrate deposit and the free gas zone forms the bottom simulating 
reflector (BSR) of acoustic waves and can be seismically detected. 

Acquiring samples with methane hydrates from sediments beneath deep water bodies is a costly procedure which 
is not practically employed for large scale surveys [11 j. Instead, the presence of hydrate is typically inferred from 
seismic data (e.g., IU2I0 or on-site salinity measurements in boreholes (e.g., [7]). The seismic data are (i) the presence 
of the "bottom-simulating reflector" which appears when hydrate dep osit touches the bottom boundary of the hydrate 
stability zone and, therefore, is underlaid by a free gas horizon 111 (Fig-ID^ an d (ii) sound speed increase in 

sediments with hydrates owing to the sediments cementation by hydrate in pores 01 211 . Both seismic techniques have 
significant limitations. For instance, the bottom simulating reflector appears only when hydrate deposit reaches the 
bottom edge of the hydrate stability zone, and the velocity increase owing to the sediment cementation does not allow 
accurate estimation of the amount of hydrate. Hence, on-site salinity measurements become an important source of 
information. 

To date, the mathematical models reconstructing hydrate deposit parameters by means of fitting the measured 
salinity profiles disregard non-Fickian contributions to the diffusion flux of salt in sediments. In the present paper we 
(i) derive relations between the profiles of the hydrate volumetric fraction in pores and the measured salinity and (ii) 
demonstrate the non-Fickian contributions to be important. 

The paper is organized as follows. In Sec. 2 we describe the transport processes in carbon-rich marine sediments 
and derive the relation between the salinity (chlorinity) and hydrate profiles. In Sec. 3 this relation is employed for 
reconstruction of the hydrate profile for a real natural hydrate deposit in the Blake Ridge hydrate province. In the 
concluding section we discuss importance of the non-Fickian diffusion and implementation of our reconstruction 
procedure. 

2. MARINE SEDIMENTS HOSTING HYDRATE DEPOSITS: TRANSPORT PROCESSES 

Real geological systems are much more uniform along two directions (say horizontal) than along the third direction 
(say vertical). Hence, we consider a one-dimensional problem with vertical spatial coordinate z (Fig.[TJ. On the field 
scale, such systems are featured by the temperature growth with depth 

T = T sf + Gz, (1) 

where T S f is the temperature of the water-sediment interface (or seafloor) and G is the geothermal gradient. 
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Methane is generated from the sediments by anaerobic bacteria. If the temperature is low enough and the pressure 
is high enough, methane forms hydrate. However, the critical pressure for the thermodynamic stability of hydrate 
depends on temperature nearly exponentially, and the hydrostatic pressure, which grows linearly with depth, cannot 
compensate the linear growth of temperature T = T s f + Gz below a certain depth zbsr- This depth zbsr is the bottom 
boundary of the hydrate deposit: below this depth, hydrate is dissociated into water and methane-gas bubbles (Fig.[TJ. 

The major part of natural hydrates of hydrocarbons on the Earth is a structure I clathrate of methane (> 99%). The 
elementary cell of an "ideal" structure I clathrate is formed by 8 molecules of CH4 and 46 molecules of H2O, i.e., the 
mass fraction of water in the clathrate Kn 2 o = 23 • 18/(23 • 18 + 4 • 16) ~ 0.866. For real hydrates the saturation of the 
structure with methane molecules is slightly less than 100% — some clathrate cages are not occupied by the methane 
molecule. In geological systems, the occupancy does not necessarily correspond to the thermodynamic equilibrium, 
because the hydrate was initially formed under the thermodynamic conditions of the place from where the hydrate was 
transported due to diverse geological processes. The relaxation rate of the occupancy is determined by the molecular 
diffusion of the methane molecules in hydrate clathrate (a solid matter) and is commensurable with rates of geological 
transport processes (100 - 1000 kyears). On the one hand, the actual local occupancy cannot be exactly evaluated 
from thermodynamical principles and the equilibrium condition, while, on the other hand, it is always higher than 



95% [ 13], i.e. very close to the ideal value, 100%. Hence, we assume the "ideal" structure of hydrate, Ku 2 o « 0.866. 

Hydrate forming in pores consumes water from the brine while salts remains in the brine. Therefore, salt con- 
centration increases and diffusion drives redistribution of the salt mass. This process determines the formation of the 
salinity (chlorinity) profile. 

Additionally, one should distinguish the on-site chlorinity, co s , and the measured chlorinity, u> s , because the drilling 
procedure results in dissociation of hydrate and the release of the hydrate water into the brine at the measurement 
site lfl~lll . Given the volumetric fraction of hydrate in pores is h(z), the mass of NaCl in the unit volume of pores 
before the dissociation of hydrate, a> s p w (l - h), equals the mass after dissociation, &> s (p w {l - h) + Kn 2 oPhh), where 
p w = 1000 kg/m 3 is the water density, pi, = 930 kg/m 3 is the hydrate density. Hence, 



v ( 1 v _ \2 -\ 

l _K M PH h + Q 



Pw 



1- H2 ° Pfl h\ = 0) s (l - kh) , (2) 

Pw 



where 



^^2^, 0.805. 



Pw 

the correction 0[(l —k) 2 h 2 ] = (9[(0.03 -h)-h] for real systems, where h rarely exceeds 7% fllH, is less than 0.002 -h 
and can be neglected. 

Since the salt is transported with pore water, we have to describe the water mass transport in the system. The 
water mass is transported as a part of hydrate, with sediments, and with the brine, by the pore water flux. The 
downward transport of sediments is significantly affected by the sediments compaction with depth |7lllJJ]; the porosity 
</> significantly decreases with depth according to the empiric law 

0(z) = o exp(- Z /L), (3) 

where L is the depth of e-folding of porosity. The mass conservation law for sediments yields 

Vj(z) = ~7~ ~77~\ v s0, (4) 

1 - 4>(z) 

where v s (z) is the sediment motion velocity. This relation is nearly not affected by the conversion of part of sediments 
into methane. Indeed, the defect of the solid matrix volume owing to methane generation is approximately AV S = 
AmcH 4 /Pi (fis ~ 2650 kg/m 3 is the sediment material density ||7|]), whereas the production of hydrate AVi, from this 
mass of methane is AmcH 4 /^HiOP/i- The ratio AVJAVf, - K^oPh/Ps ~ 0.04 is small and thus AV S related to methane 
generation is negligible. 

For water, the mass conservation law reads 

d d d 

—(0(1 - h - b)p w + 4>hK H , ph) = —^-ipwiiw) - -^-(<t>hK H ,op h v s ) , (5) 
ot az oz 
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where b is the volumetric fraction of bubbles in the pore volume (h — beyond the hydrate stability zone, b — 
within it), u w is the brine filtration velocity. Following Davie and Buffett D7l|8[], we consider a steady-state situation 
and set time derivatives to zero. Hence, 

p w u w (z) + cf>(z)h(z)K H2 op h v s (z) = p w u w (0) 

(hydrate is not present close to the water-sediment interface, h(Q) = 0), and, substituting v s (z) from Eq. ©, we find 

u w = u w0 - k(f>h- — ^j-Vsq . (6) 
1 -0 

Given the brine filtration velocity u w is known (Eq. (0), one can evaluate the salt transport from the mass conser- 
vation law; 

d d d 

—(0(1 - h - b)p w u) s ) = - — (p w co s u w ) - — 7. 5 ,difr • (7) 
ot dz az 

The diffusive flux J s ,m of the NaCl mass is contributed by the Fickian molecular diffusion flux and non-Fickian 

diffusion fluxes — the thermal diffusion and gravitational segregation (importance of which for the gas transport in 

geological systems under consideration was previously demonstrated [5]). The diffusive flux reads ifTil BIl 

/ s ,diff = -#0(1 -h- b)D s p w u s ^ In cj s + a s -^- In T - ^ j . (8) 

The following notations are introduced here: 

• X is the tortuosity factor, which characterize the effect of the pore morphology on the effective diffusivity of species 
in pore fluid. For our system x — 0-75 II 1511 . 

• D s is the molecular diffusion coefficient in bulk brine. 

• a s is the thermodiffusion constant. 

• g = 9.8 m/s 2 is the gravity. 

• The universal gas constant R = 8.3 14 J/ (mol K). 

• fi s = //NaCl - Np.H 2 o is the effective molar mass of the pair of ions Na + and CF in the aqueous solution, is 
the number of water molecules in the volume occupied by this pair in the solvent, which can be evaluated from the 
dependence of the solution density on its concentration (see, e.g., Appendix A in ||5|]): 



A<NaCl <9p so lution 

Ps = 



42 g/mol . 



Pw da> s 

In d] the thermodiffusion constant a s and the molecular diffusion coefficient D s were measured for a seminormal 
aqueous solution of NaCl. For the thermodiffusion constant a sign inversion was observed near T{ = 12°C. In the 
temperature range typical for our system, T e (275 K, 305 K), the temperature dependencies of a s and D s are strong 
and well represented by expressions 

a s * 0.0246 K'^T-Ti) and D s * 6.1 • 10~ 10 exp[0.0371 Kr'(r - 273.15 K)] m 2 /s . 

Hydrostatic pressure in marine sediments, which is up to several hundreds atmospheres, is not strong enough to affect 
the diffusion constant of chemicals in water. 
For a steady state, Eqs. (0 and © yield 

co s u w -#0(1 -h - b)D s u) s |^ ln^ + a s ^ In T - ^pj = o).„u m + J sm „ , (9) 

where (j s , = co s (z*), u w * = u w (z*) = u„o, J s ,m* = Js.diffiz*), and z* is a certain depth deep below the hydrate stability 
zone (see Fig.Q}. Hereafter, the sign "*" indicates the value at depth z«. 

Since h <k 1, we restrict our consideration to the linear in h (and b) approximation. Substituting Eqs. (f2]i and © 
into Eq. (|9|i we can find 

dh 

— +y( z )h = f(z), (10) 

dz 
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water-sediment interface temperature, Eq. ([T] 


2°C 


Refs. [11, 7J 


G 


geothermal gradient, Eq. © 


35°C/km 


Refs. \VIT\ 


<f>(0) 


seafloor porosity, Eq. Q 


0.69 


Ref. CJ 


L 


e-folding depth of porosity, Eq. (0 


2 km 


Ref. [7J 



where 

1/1-00 «wo\ . 1/1 <3&,s G\ <t>*D st I 1 <9a>.« G*\ m,,o /, S) s 
r(z) = 77T h r v ^-— , /(z) = - 7 +7T7T ~ 5~ + £*7f- + 



^ D s \ 1 - ' cp j k\<l> s dz Tj kcf>D s \<jj st dz T„j kx<pD s \ <1> S j ' 

and the parameter 

P = aAT) -RG 

characterizes the strength of the non-Fickian flux. 

With Eq. (flOl i. one can reconstruct the hydrate profile h(z) from the measured chlorinity profile ib s (z). Although 
one can write down an analytical solution to the problem (flOt 



h(z) - C /(:,)<■ ! y> : " L < lz ] . 
Jo 



numerical integration of Eq. (flOt is more convenient for data analysis in practice. Remarkably, the relation between 
h{z) and Cj s {z) does not involve quantitative data on the process and history of the formation of hydrate deposit and 
the process of generation of methane from sediments. 



3. SALINITY PROFILE ANALYSIS AND HYDRATE PROFILE RECONSTRUCTION 

We demonstrate application of our results to the analysis of one of the most important marine hydrate provinces — 
the Blake Ridge. For the Ocean Drilling Program site 997, on the Blake Ridge, extensive data have been acquired, 
including hydrate samples and the on-site salinity (chlorinity) measurements [11] (Fig.|3a)). The reported parame- 
ters for this site are presented in Table [1] We have two parameters which are not imposed by the results of direct 
measurements: sedimentation rate v,o and filtration velocity u w q. 

In natural systems, the hydrate deposit cannot be in touch with the water-sediment interface, because aqueous 
methane concentration in the water body above the sediments is zero and hydrate must dissociate. Moreover, in seas, 
methane is oxidized by sulfates, which are present in sea water, and its concentration is zero within the so-called sulfate 
reduction zone, which typically expands approximately 20 m below the water-sediment interface [§[]. Hence, hydrate 
should not be present in a quiet extended upper part of the hydrate stability zone. With Eq. ([Tol l, the absence of hydrate, 
h = 0, requires f(z) - next to z = 0. The function f(z) is independent of VsO and we can set it to zero for small z by 
tuning u w q. With the chlorinity profile plotted in Fig.|2ja), this procedure yields u w o = -(8 + 0.5) cm/kyear (the flux 
is negative, i.e. ascending). One can see that for u w o = -8 cm/kyear, h — down to depths slightly over 100 m, while 
for different filtration velocity it deviates from zero next to z — m. Furthermore, with fixed u w o, the reconstructed 
amount of hydrate below the hydrate stability zone (for z > 450 m in Fig.|2j depends on v s o monotonically ; it vanishes 
for Vjo = (9 ± 0.5) cm/kyear (Fig.|2b)). The hydrate profile plotted in Fig. 0b) with the red solid line is the final result 
of the reconstruction of the hydrate profile from the measured chlorinity profile plotted in Fig.|2ja). 

It is noteworthy that our reconstruction procedure is free of uncertainties in parameters: all but two parameters are 
available from direct measurements and these two parameters are strictly imposed by two inevitable inherent features 
of the hydrate profile. 
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Figure 2: (Color online) (a): Measured chlorinity profile (solid circles) and smoothed chlorinity profile used for calculation of the hydrate profile 
(red solid line) for the site 997 of the Ocean Drilling Program Hill , (b): Hydrate profile reconstructed from the chlorinity data with Eq. i\0\ is 
plotted with the red solid line (parameters are specified in Tableland n„,o = -8 cm/kyear, v,o = 9 cm/kyear). For demonstration, we plot a formal 
hydrate profile for purely Fickian diffusion flux of the same strength (blue dashed line). Non-Fickian contributions are obviously non-negligible. 
Indeed, the latter profile significantly deviates from the former one and possesses unphysical features: negative values of the hydrate volumetric 
fraction and non-zero (negative) amount of hydrate beyond the zone of the thermodynamic stability of hydrate. 



4. CONCLUSION AND DISCUSSION 

In this paper the transport of water and salt have been considered for marine sediments hosting natural hydrate 
deposits. The mathematical description employed accounts for 

• non-Fickian diffusion of NaCl, and 

• temperature dependence of the molecular diffusivity. 

We have demonstrated the crucial importance of the both, whereas they are disregarded in the literature on the mod- 
elling of marine hydrate deposits (e.g., Sill Got). Based on this consideration, we have derived the relation between 
the measured salinity (chlorinity) profile and the hydrate profile. Application of this relation has been demonstrated 
for a real hydrate deposit (Fig. [2). 

The solution, we found for this "reverse engineering" problem, does not involve quantitative data on the process 
and history of the formation of hydrate deposit and the process of generation of methane from sediments. This is 
an important feature of our results because previously, in the literature, closed models of hydrate formation involve 
particular assumptions on the generation process (e.g., (0, [Hi)- In these studies the entire model is tested against 
the measured salinity profile, while we can see that only the current hydrate profile determines the salinity profile. 
Moreover, the sedimentation rate is unambiguously imposed by features of one of these profiles, whereas it has been 
previously indirectly inferred from geological data. 
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Importantly, our reconstruction procedure is free of uncertainties in model parameters: all but two parameters — 
sedimentation rate v s q and filtration velocity u w q — are available from direct measurements. These two parameters are 
strictly imposed by two inherent features of the hydrate profile: the absence of hydrate (i) close to the water-sediment 
interface and (ii) beneath the hydrate stability zone. 
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